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Abstract —Fundamental theory on battery-powered cyber¬ 
physical systems (CPS) calls for dynamic models that are able 
to describe and predict the status of processors and batteries at 
any given time. We believe that the idealized system of single 
processor powered by single battery (SPSB) can be viewed as 
a generic case for the modeling effort. This paper Introduces a 
dynamic model for multiple aperiodic tasks on a SPSB system 
under a scheduling algorithm that resembles the rate monotonic 
scheduling (RMS) within finite time windows. The model contains 
two major modules. The first module is an online battery capac¬ 
ity model based on the Rakhmatov-Vrudhula-Wallach (RVW) 
model. This module provides predictions of remaining battery 
capacity based on the knowledge of the battery discharging 
current. The second module is a dynamical scheduling model 
that can predict the scheduled behavior of tasks within any finite 
time window, without the need to store all past information about 
each task before the starting time of the finite time window. 
The module provides a complete analytical description of the 
relationship among tasks and it delineates all possible modes of 
the processor utilization as square-wave functions of time. The 
two modules i.e. the scheduling model and the battery model are 
integrated to obtain a hybrid scheduling model that describes 
the dynamic behaviors of the SPSB system. Our effort may 
have demonstrated that through dynamic modeling, different 
components of CPS may be integrated under a unified theoretical 
framework centered around hybrid systems theory. 

I. Introduction 

Modern society is characterized by the pervasive applica¬ 
tions of embedded computing systems powered by batteries. 
Use of batteries endows computing systems with mobility that 
is appreciated by consumers, industries and governments. Cell 
phones, portable music/video players, and unmanned airplanes 
flying over a battle zone are among the most noticeable 
applications. Computing systems supported by batteries are 
usually designed following specialized guidelines to reduce 
power consumption, achieving lower maintenance and longer 
operation time. 

Despite of the great achievements in battery-powered com¬ 
puting, better theoretical understandings of the interactions 
between computing systems and batteries will be appreciated 
by applications where power is very limited or battery life is a 
dominant constraint on systems design. One such application 
is in underwater robotics, where missions may last more than 
a year. Although appear to be simple, batteries are complex 
physical-chemical systems. Hence fundamental research in 
battery powered embedded computing systems may be viewed 


as part of the evolving cyber-physical systems (CPS) theory 
(c.f. recent publications HI, ||2l, E], H, jS), IQ) that aims 
to develop novel foundations to understand complex physical 
systems controlled by modern computers. 

The complexity of battery discharging behaviors is noticed 
in literature e.g. Q, 0. Battery modeling aims to simulate 
these behaviors by computational models 0, ifTOl . To support 
theoretical cyber-physical systems design, a class of analytical 
models are desired. This is because comparing to physical 
Ga, empirical IH, and circuit based lEi, ca models, 
analytical models are theoretically tractable while providing 
sufficient accuracy. 

The interactions between batteries and computing systems 
have been studied by some researchers. Power consumption 
of computing devices is proportional to lfl5l . hence can 
be reduced by dynamic voltage scaling (DVS). Yao et al. M 
propose one of the first DVS-aware scheduling algorithms. 
Rowe et. al. iflTll proposed a scheduling method to combine 
idle period of tasks so that the processor can be put to sleep 
to save power. Some recent results on power management in 
sensor network applications are presented by Ren et. al. m. 
In these works the behaviors of batteries are simplified as 
ideal voltage sources whose life only depends on the average 
discharge current. The approach taken by Jiong et. al. Cl, 
ll20ll uses an empirical battery model and assumes known 
battery discharge profiles for computing tasks. 

We believe that a dynamic model needs to be established 
that is able to predict battery capacity based on the discharge 
current determined by a feedback control law or feedback 
scheduling algorithm. The magnitude and pulse width of 
such discharge current can not be determined beforehand. 
Therefore, in previous work mi, we introduce a dynamic 
battery model that describes the variations of the capacity of 
a battery under time varying discharge current. This model is 
input-output equivalent to the Rakhmatov-Vrudhula-Wallach 
(RVW) model proposed by Rakhmatov et. al. m which 
has been shown to agree with experimental results and has 
demonstrated high accuracy in battery capacity prediction and 
battery life estimation. Besides the accuracy, the RVW model 
is also simple to use. It can characterize any general battery 
with only two parameters a and /3, which can be estimated 
by a experimental method introduced in ll22l . However, the 
application of the RVW model requires that the discharge 
current, as a function of time, is known. Our dynamic model 
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allows battery capacity prediction for feedback control laws 
and online scheduling algorithms. 

This paper is inspired by the need to completely describe 
the possible battery discharge profiles for a processor that is 
running a finite number of aperiodic real-time tasks. Knowing 
such a discharge profile allows us to predict remaining battery 
capacity at any given moment without seeking the help from 
simulation programs. We found that the problem can only be 
solved by establishing a dynamic model that integrates both 
battery models and task scheduling models. However, we are 
unable to find existing models from the literature that fits in 
such requirements. 

Therefore, we start our modeling effort by considering a 
single processor single battery (SPSB) model. This model 
may appear naive. However, we find it is very challenging 
to establish a complete analytical description for the behav¬ 
iors of multiple tasks. Readers in doubt about this difficulty 
may consider giving an explicit mathematical formula for 
the processor utilization plots generated by real-time systems 
simulation tools such as the TrueTime 1^ . Such difficulty 
has been mostly avoided by the works on schedulability 12^ . 
Il25l . 12^ because only the worst case scenario, which are the 
critical time instances, are considered. 

With the help from hybrid systems theory evolving in the 
control community Il27l . ED, we obtain a class of hybrid 
dynamical systems models that are able to describe the battery 
discharge profile when multiple independent aperiodic tasks 
are scheduled under a fixed priority scheduling algorithm on 
a single processor. When integrated with the dynamic battery 
model, theoretical predictions of battery capacity for an SPSB 
system can be produced online at any given time. Comparing 
to the literature reviewed, we believe our contribution is novel 
and may provide foundation for developing novel battery 
aware scheduling for cyber-physical systems. 

The paper is organized as follows. Section Hill describes the 
SPSB system we investigate. Section HJ introduces a dynamic 
battery model that is used to determine the remaining battery 
capacity. Section m derives the dynamic models for two 
aperiodic tasks on a single processor. The dynamic model 
also enables us to find analytic descriptions of the processor 
utilization waveforms in Section |V] Section |Vl] extend the 
conclusion of two tasks to multiple tasks. The models are 
verified and demonstrated by simulations in Section IVIII 
Conclusions are provided in Section IVIIII 

H. A Dynamic Battery Model 

In this section, we briefly review the dynamic battery model 
established in our recent work m as a workshop paper. We 
model the discharge dynamics within one discharge cycle of 
(possibly rechargeable) batteries. The discharge cycle starts 
when a battery is fully charged and ends when the battery can 
not support the demand for discharge current; no recharging 
is allowed during the cycle. According to battery and VLSI 
design literature e.g. Q, 0, the discharge current is supported 
by the change of concentration of electrolytes near the elec¬ 
trodes of a battery. When the concentration at the electrodes 
drops below a certain threshold, a battery fails to support 


discharge current, causing failure to devices it supports. At 
this moment, the discharge cycle has to stop and the battery 
needs to be recharged or replaced. Note that there may be a 
significant amount of electrolytes left when a discharge cycle 
ends. When the battery is discharged under a pulsed discharge 
current, during the idle time when current is interrupted, the 
diffusion process increases electrolyte concentration at the 
electrodes. This produces the recovery effect that makes the 
battery appears to have regained portions of its capacity. 

The Rakhmatov-Vrudhula-Wallach (RVW) model in ll22ll is 
derived from solving the diffusion equations governing the 
electrolytes motion within a battery. The model captures the 
recovery effect effectively. As mentioned in the introduction, 
in order to use the RVW model for CPS design, we derived 
a dynamic model that is input-output equivalent to the RVW 
model. The key improvement is to introduce a set of state 
variables, denoted by x = [xo,a;i, where m is chosen 

according to accuracy requirement. In most cases, m < 10 is 
accurate enough. Then the dynamic model is given by a linear 
system: 

X = Ax-|-b2(f) 

y = c'x. (1) 

Here, the matrix A is a diagonal matrix i.e. A = 
diagjO, —Ai,..., —Am} where Xj = (3 * ff for j = 1, 2,..., m. 
The input i{t) > 0 is the discharge current. The vector 
b = A [1, 2, 2,..., 2]' and the vector c = [1,1,1,...,!]'. 
As mentioned in the introduction, a and /3 characterize a 
battery and they can be estimated by the experimental method 
introduced in ll22ll . The output variable y measures the total 
capacity loss for a battery. It contains two parts, the permanent 
capacity loss represented by xq, and the temporary capacity 
loss measured by xi,X 2 , At the beginning of the 

discharge cycle the initial condition for x is the zero vector, 
hence y = 0. When y = 1, the discharge cycle ends. The first 
state variable xq is the integration of the discharge current 
i(t), hence is the effective discharge delivered to the circuits 
outside the battery. One can see that the actual discharge y 
contains more than xq. The states xi, X 2 ,..., Xm are temporary 
capacity losses that are always positive since i{t) > 0 . Under 
an impulsive current, when iff) = 0 , Xj(l < j < rn) 
decreases because —\j < 0 . This captures the recovery effect. 
More details about this model and results regarding optimal 
discharge profiles are presented in our work ED. 

III. Scheduled Behaviors of Aperiodic Tasks 

In this paper, we consider the case where aperiodic tasks 
with hard deadlines scheduled on a SPSB system. Tasks are 
scheduled using the following scheduling algorithm: the tasks 
with shorter request intervals are assigned higher priorities, 
and the higher priority tasks can preempt the lower priority 
tasks. We are interested in the scheduled behaviors of 
Ti,..., Tat within a finite time window, which starts from fo 
with length L. to can be any time instant. Since the request 
intervals of aperiodic tasks are not constant, the priority 
assignments among aperiodic tasks may change during the 
finite time window. To avoid such cases, we enforce the 
following assumptions on the length of time window L. 
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Assumption 3.1: Within [to, to + L], the maximum request 
interval of is smaller than the minimal request interval of 
Ti+i, where i is selected such that the priority of is always 
higher than the priority of r^+i. 

When to changes, the length of time window and the priority 
assignments may change as well. 

Assumption 3.2: We assume that for any time instant t, 
there exists to and L such that Assumption 13. II is satished. 

Assumption 13.11 guarantees coherence in priority assign¬ 
ments during the finite time window starting from to. As¬ 
sumption 13.21 guarantees that there is no gap between the 
proper time windows. It is of course trivial to see that periodic 
tasks satisfy these assumptions within + L]. In fact, the 
scheduling algorithm becomes the Rate Monotonic Scheduling 
(RMS) when the time window L goes to inhnity. Our main 
goal here is to model the behaviors of the tasks, not proposing 
a new scheduling algorithm. 

Since such an modeling effort has not been performed 
in previous literature, we need to introduce new notations 
that depict the task behaviors. Suppose Tq, is an independent 
aperiodic task. We use ta{n) to denote the request time of 
the n-th instance of Tq,; Ca{n) to denote the computing time 
of the n-th instance of and Ta{n) to denote the request 
interval of the n-th instance of r„, which is dehned to be the 
interval between the request time of the n-th instance of Tq, 
and the request time of the (n + l)-th instance of r^, i.e. we 
have Ta{n) = ta{n -f 1) — ta{n). 

Regardless of the choice of to, the scheduled behaviors of 
ri,...,rjv within [to, to -f L] is affected by both the tasks 
requested within + L] and the tasks requested before 

to. To predict the scheduled behaviors within + T] 

by real-time simulation tools such as Truetime, we have to 
know all task information within [0, to + L] and simulate from 
time 0. However, keeping all task information is costly in the 
real application and simulating from the beginning is time 
inefficient. We discover that this problem can be solved by 
using a dynamical model, instead of storing a large amount of 
data, we just need to update the value of some state variables. 
Furthermore, the dynamic model is able to determine the state 
of the processor at any given time within [to, to + L]. 

We develop a recursive algorithm to determine the state 
of the processor for multiple aperiodic tasks. Our goal is to 
hnd out when the processor is occupied within + L]. 

The result is a function of time $ : t {0,1} where for 
t G [to, to + L], $(f) = 0 if the processor is free, and <i>(f) = 1 
if the processor is busy. 

Furthermore, we assume that when $(f) = 1, a constant 
current will be drawn from the battery and when <l>(f) = 0, 
the current vanishes. For simplicity, we ignore the possible 
power difference among tasks and the transients in the battery 
discharge current. If $(f) is known, we are able to describe 
the current i{t) drawn from the battery by the processor. And 
then, using the dynamic battery model, we are able to predict 
the remaining battery capacity. 

IV. A Dynamic Scheduling Model eor Two Tasks 

In this section, we establish a dynamic model for the 
scheduled behaviors of two independent aperiodic tasks on a 


single processor within [fo,fo + L]. Later in this paper, we 
will develop a recursive algorithm to extend the model to 
multiple tasks. Consider two tasks and Tp, where a = 1 and 
(3 = 2, scheduled under the scheduling algorithm introduced 
in Section [111] Suppose Assumptions 13.11 and 13.21 are satisfied, 
we have the following claim; 

Claim 1: Tq, is assigned higher priority than Tp within 
[fo: to + L]. 

There is at most one instance of Tp requested within each 
request interval of r^. 

We determine a set of ’’state variables” that completely 
determine the scheduled behavior of and Tp at any given 
time within [to, to -I- L]. Then transition rules between these 
states from the current time to all future times are described 
by a set of nonlinear difference equations.The work presented 
in this section centered around these two themes—states and 
transition rules. 

A. Exploring Task Relationship 

In this subsection, we model the dynamics for the request 
time of Tq and Tp, and explore task relationship between Tq 
and Tp. 

We hrst select the request time of Tq and Tp as two state 
variables. Within the hnite time window, ta{n) and tp{m) 
satishes the following equations 

ta{n + l) = ta{n)+Ta{n) 
tp{m -I- 1) = tp{m) + Tp{m) 

where n and m are the index numbers for the instances of Tq 
and Tp requested within the finite time window. For example, 
if Tq has a fixed period, n ranges from 1 to [L/Ti\ where Ti 
is the period for Tq; and if Tp has a fixed period T2, m ranges 
from 1 to [L/T2J. 

Next, we explore task relationship between Tq and Tp by 
introducing another two state variables y(n) and z{n) . y{n) 
helps describe the request time of the hrst instance of Tp 
requested after the n-th instance of Tq, which is denoted by 
tp{y{n)). The value of y{n) indicates the index numbers for 
the instances of Tp requested within the hnite time window. If 
Tp has hxed period T2, then y{n) ranges from 1 to [L/T2J. 
z{n) describes the phase difference between the request time 
of the n-th instance of Tq, and the request time of the ?/(n)-th 
instance of Tp, i.e. 

z{n) =tp{y{n))-ta{n). (3) 

To determine the transition rules for y(n) and z{n) from n 
to n -F 1, we need to consider two cases. 

Case 1 This case happens when the phase difference 
satishes 

0 < z{n) < Ta{n) (4) 

In this case, the y{n)-t\\ instance of Tp is requested within 
\ta{n),ta{n + 1)]. Thus, the y{n)-t\\ instance of Tp is the last 
instance of Tp requested before ta{n + l). On the other hand, 
since y{n + l) is dehned to be the index number for the hrst 
instance of Tp requested after ta{n + l), the {y{n + l) — l}-th 
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instance of rp is also the last instance of Tp requested before 
ta{n + 1). Therefore, we have 

y{n + l) - 1 = y{n) (5) 

In this case, z(n) will update in the following ways 

z{n + l) = t/ 3 {y{n + 1)) - ta{n + 1) 

= t/ 3 {y{n) + 1 ) - ta{n + 1 ) 

= z{n)+T 0 {y{n))-Ta{n) 

where (|2]i, 0 and Q have been applied. This implies that 
the phase difference between the two tasks has increased by 

Tp{y{n)) -Ta{n). 

Case 2 This case happens when the phase difference 
satisfies 

z{n) > Tain) 

In this case, the y{n)-th instance of is requested after 
tain + 1). Thus, the y{n)-th instance of is the first instance 
of Tp requested after tain+1). On the other hand, the ?/(n+l)- 
th instance of Tp is also defined to be the first instance of rp 
requested after tain + 1). Therefore, we have 

yin + 1) = yin) (6) 

In this case, z(n) will update in the following ways 

zin+l) = tpiyin + 1 ))-tain + 1 ) 

= tpiyin))-tain + 1) 

= zin)-Tain). 


to compute the last instance of rp requested before tain), i.e. 
the {yin) — l}-th instance of rp, . 

To determine the value of i?(n), we need to consider two 
possibilities. 

1) If 

0 < Pin) < Tain)-Cain), (9) 

which implies that the execution of the [yin) — l}-th 
instance of rp is finished within [tain),tain + l)], then 
i?(n) = Pin). 

2) If 

Pin) > Tain) - Cain), (10) 

which implies that the execution of the {yin) — 1}- 
th instance of rp cannot be finished before tain + 1), 
then the idle time of processor within the interval of 
\tain), tain + 1)] will be allocated to compute rp. Thus, 
Pin) = Tain) - Cain). 

As a summary, we have 

i?(n) = min{T(j(n) — Cain), Pin)}. (11) 

We can see that i?(n) solely depends on Pin), but i?(n) 
will be more convenient to use to derive processor utilization 
waveforms in Section |V] 

Next, we determine the transition rules for Pin). There are 
two cases to consider. 

Case 1; the phase variable z(n) satisfies 

0 < zin) < Tain) (12) 


It can be seen that the phase difference decreases by Tain). 

In summary, y(n) and z(n) are described by the following 
equations: 


yin + 1 ) 


yin) + l if 0<zin)<Tain) 
yin) if zin) > Tain) 


( zin)+ Tpiyin))-Tain) 

I {if 0 < zin) < Tain)} 

I zin)-Tain) 

[ (if -zin) > Tain)} 


which implies that the j/(n)-th instance of rp is requested 
within [tain), tain + 1)], as illustrated by Fig. [T] 




B. Modeling the residue 

The state of processor after a given time point will be 
affected by the instances of Tq and rp requested before the 
given time point. To track the status of Tq, and rp, we introduce 
a term called the residue. The residue P(n) is defined to 
be the amount of time that is needed after the time instant 
tain) to finish computing the last instance of rp requested 
before tain). More specifically, since the j/(n)-th instance of 
Tp is defined to be the first instance of rp requested after 
tain), the {yin) — l}-th instance of rp is the last instance 
of Tp requested before tain). As a result, P(n) describes the 
remaining time needed to finish computing the {yin) — l}-th 
instance of Tp after tain). This P(n) will be another state 
variable, in addition to tain), tpin), yin) and zin). We also 
introduce an auxiliary variable P(n) that describes how much 
time within the interval of [tain), tain + 1)] can be allocated 


Fig. 1. Case 1. The y ( n )- th . instance of is requested within [fa (n), fa in -\- 
1)] 

In Fig. [T] we can see that if 0 < zin) < Cain), then 
the execution of the ?/(n)-th instance of Tp will start from 
ta in) + Ca in) since the execution of the n-th instance of Tq 
will postpone the ?/(n)-th instance of Tp. If Cain) < zin) < 
Tain), then the execution of the 2 /(n)-th instance of Tp will 
start from tain) + zin). Thus, the execution of the ?/(n)-th 
instance of Tp will start from 

max{fQ,(n) + Cain),tain) + zin)}. 

The time allocated for the execution of the y(n)-th instance 
of Tp within [tain), tain + 1)] is, 

tain + 1) — max{fa(n) + Cain), tain) + zin)} 

= Tain) — max{Cain), zin)} ( 13 ) 
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Therefore, the residue for the ?/(n)-th instance of Tp to be 
computed after ta{n + 1) is 

Cpiyin)) - (Ta{n) - max{C'a(n), 2 ;(n)}). (14) 

The value of P(n + 1) can not be negative. Therefore 

P(n+1) = max{ 0 ,C', 3 (j/(«))-(T'a(?T-)-max{C'a(n),z(n)})}. 

(15) 

Case 2; When the phase variable z{n) satisfies 


We always assume that tc{l) = to- If fo 7^ we will 
model a new instance of as the first instance of requested 
within [fo, fo + L]- The new instance has the following task 
characteristics: 

tM)=to CM) = Po. T„(l) =f,(l)-fo. (19) 

According to the definition of tp, P, y and z, we have 
tpil)=tp y(l) = l z{l) = tpil)-tM) P{^)=Pp 


z(n) > Ta(n) (16) 

which implies that the y(n)-th instance of is requested after 
ta(n + l), the {y(n) — l}-th instance of Tj 3 is requested before 
ta(n), and no instance of Tp is requested within [ta(n),ta(n+ 
1)], as illustrated by Fig. |2 . 



C„(n) 




k(n+l) 


z(n) 


z(n)> (n) 

Fig. 2. Case 2. The y(n)-th instance of is requested after ta(n + 1) 


In this case, the {y(n) — l}-th instance of rp is the last 
instance of rp requested before ta(n). By definition, we know 
that P{n) describes the remaining time to compute the {y{n) — 
l}-th instance of rp after ta{n). In addition, the {y{n) — l}-th 
instance of rp is also the last instance of requested before 
ta{n + 1) and thus P{n + 1) describes the remaining time to 
compute the {y{n) — l}-th instance of Tp after ta{n + 1). 

From Fig. 121 we can see that the time that can be allocated to 
compute the —l}-th instance of Tp within [ta{n),ta{n+ 
1)] is Ta{n) — Ca{n). Therefore, the remaining time for the 
{y{n) — l}-th instance of rp to be computed after ta{n) is 

Pin) - (Tain) - Cain)). 

Thus the value of P(n + 1) is given by the following: 


Pin + 1) = max{0, P(n) - iTain) - Cain))}. (17) 


According to O and (O, we can draw a conclusion on 
the transition rules for the state variable Pin) from n to n+1. 


Pin + 1 ) = < 


ma,x{0,Cpiyin)) - iTain) 

-max{C'a(n),z(n)})} 
if 0 < zin) < Tain) 
max{0, Pin) — iTain) — Cain))} 
if zin) > Tain) 

(18) 


C. Initialization of State Variables 

In this subsection, we will discuss the initialization of the 
state variables within [to, to + L]. 

Suppose ta and is the request time of the first instance 
of Ta and Tp requested within [fp, to + L] , Pa and Pp is the 
residue for the last instance of and Tp requested before to 
to be computed after to- Pa and Pp need to be retrieved from 
the memory of the system. 


As we can see, the residue P serves to connect the state of 
processor before to with the state of processor within [fo, to + 
L]. Our model has the advantage over the classical scheduling 
analysis because instead of storing for all task behaviors made 
before to, we only need to record the residue of each task and 
to start the model. Therefore, we can analysis the scheduled 
behaviors of tasks starting from any time and within any finite 
time window as long as Assumptions 13.11 and 13.21 are satisfied. 

V. Processor Utilization Waveform of Two Tasks 

Suppose Ta and Tp are schedulable. With the states and 
transition rules for the scheduling algorithm determined, we 
compute a function $(<) that describes the states of the pro¬ 
cessor as either free ($(f) = 0) or busy ($(f) = 1). To achieve 
this goal, we study the states of processor within each request 
interval of Ta, i.e. [tain), tain + 1)]. The result is a function 
of time t —s- {0,1} where for t S [tain),tain + 1)], 

T’r!,(i) = 0 if the processor is free, and $„(f) = 1 if the 
processor is busy. If is known for all n, we are able 

to describe $(f) within + L]- Note that we ignore the 

scheduling transients that may exist for the processor, then the 
graph of such a function resembles a square wave. Although 
the space of possible square wave functions on a compact 
interval is infinite dimensional, we find that for the scheduling 
policy introduced in Section Hn] the waveforms only have two 
modes. 

Case 1 This describes the mode when the phase variable 
zin) satisfies 

0 < zin) < Tain), (20) 

which implies that the j/(n)-th instance of Tp is requested 
within [tain), tain) +Tain)]. 

There are two subcases. 

Case 1.1: If the phase variable z(n) satisfies 

0 < zin) < Cain), (21) 

which implies that the j/(n)-th instance of Tp is requested 
while the processor is computing Ta, i.e. tpiyin)) G 
[tain), tain) + Cain)], as illustrated by the left picture in 

Fig. [II 

It can be observed in Fig. [T| that 

$„(f) = 1 if t G [tain),tain) + Cain)]. (22) 

The {yin) — l}-th instance of Tp must have been finished 
before tain) as a result of schedulability requirements. Thus, 
i?(n) = 0. In order to derive a unified formula for later, 
we notice that (l22ll can be rewritten as 

= 1 if t G [tain), tain) + Cain) + Pin)]. (23) 

















6 


The execution of the y(n)-th instance of Tp will start from 
tain) + Cain). The finishing time of the j/(n)-th instance of 
Tp within [tain), tain + 1)] have two possibilities: 

1) If 

C^pivin)) < Tain) - Cain) (24) 

which implies that the execution of Tp will be finished 
before tain + 1), then 

$n(f) = 1 if 

t e [tain) + Cain), tain) + Cain) + Cpiuin))]. 

(25) 

2) If 

Cpivin)) > Tain) - Cain) (26) 

which implies that the execution of Tp will occupy the 
idle time of processor within the interval [ta in ), (n + 
1)], then 

$„(f) = 1 if 

t e [tain) + Cain), tain) + Cain) + iTain) - Cain 

ill) 

As a summary, the execution time of the 2 /(n)-th instance 
of Tp during [tain), tain + 1)] is 


As a summary, a conclusion can be drawn from (l29l l and 
(l3Tt that during [tain), tain + 1)], the processor state is 

$„(f) = l if 

t e [tain), tain) + Cain) + Rin)] IJ [<„(«) + 

max{C'ti(n), zin)}, tain) + max{C'Q-(n), zin)} 
-\-m.m{Cpiyin)),Tain) - max{C'a(n), ^(n)}}]. 

Case 2 This describes the mode when the phase variable 
zin) satisfies 

zin) > Tain). (35) 

which implies that the y(n)-th instance of Tp is requested after 
+ 1), as shown in Fig. |2] 

In Fig. |2] only the {yin) — l}-th instance of Tp and the 
n-th instance of will affect the processor state within 
[tain),tain + 1)]. In this case, the {yin) — l}-th instance 
of Tp is the last instance of Tp requested before tain). Thus, 
the execution of the {yin) — l}-th instance of Tp and the n-th 
instance of Tq, will take up the time 

Rin) + Cain) 

which implies that the processor state during [tain), tain+1)] 
is 


$„(f) = l if 

t G [tain)Cain), tain)Cain)-\- (28) 

m.m{Cpin). Tain) - Cain)}]. 

According to (|2^ and (l28T l. the processor utilization during 

[tain),tain + 1)] is 

$„(f) = 1 if 

t € [ta in), ta in) + Ca in) + Rin)] [J [ta in) + Ca in), 
tain) +Cain) +min{Cpin), Tain) - Cain)}]. 

(29) 

Case 1.2: If the phase variable zin) satisfies 

zin) > Tain) (30) 

which implies that the j/(n)-th instance of Tp is requested 
within [tain) + Cain), tain + 1)] , as shown in the right 
picture of Fig. |2] 

The arguments for the {yin) — l}th instance of Tp still 
holds. Hence 

= 1 if t&[tain),tain)+Cain) + Rin)]. (31) 

The execution of the y(n)-th instance of Tp will start from 
tain) + zin). The finishing time of the ?/(n)-th instance of 
Tp within \tain),tain + 1)] have two possibilities. Using the 
similar method for (l28T l. we can show that the execution time 
of the y(n)-th instance of Tp during [tain), tain + 1)] is 

4>„(f) = 1 if 

t & [tain) + zin), tain) + zin)+ (32) 

min{C',3(2/(n)), Tain) - zin)}]. 

According to (lUT i and (l32l i. we know that the processor 
utilization during [tain), tain + 1)] is 

4>„(f) = 1 if 

t € [tain), tain) + Cain) -I- Rin)] |J [tain) + zin), 
tain) + zin) + m.m{Cpiyin)), Tain) - zin)}]. 

( 33 ) 


= 1 if t e [tain), tain) + Cain) + Rin)] . (36) 

According to (l34l i and (l36l l. the processor utilization within 
each request interval of Tq, are depicted in Fig. |3] 


C„(n)+R(n) min{C^(y(n)),T„(n)-max{C„(n),z(n)}} 

Case1 ^ 

C„(n)<z(n)<U(n) _ I 2 | _□ 

k(n)+maxp„(n),z(n)} t^(n4-I) 


Case2 ^ 

z(n) > (n) □ 

U(n) 


C„(n)+R(n) 


F(n+1) 


Fig. 3. Processor Utilization Waveform within [ta{n), tain + 1)] 


VI. Model For Multiple Tasks 

In Section |IV] and |Vl we derive the processor utilization 
waveform of two tasks from a dynamic scheduling model. 
In this section, we develop a recursive algorithm to extend 
the modeling effort from two tasks to multiple tasks. Suppose 
Ti, ...,tm are schedulable. The processor utilization waveform 
of Ti, ...,tn can also be derived recursively in the following 
ways 

1) Initialization: In the first recursion, we use Ta to denote 
Ti and Tp to denote T 2 . Then, we derive the processor 
utilization waveform of Ta and Tp according to (l3^ and 
(l36l) . and model the waveform into a single task Tcmb- 

2) Iteration: In the n-th iteration (2 < n < — 2), we 

use Ta to denote the Tcmb derived from the (n — 1)- 
th iteration and Tp to denote Then, we model the 

processor utilization waveform of and Tp into a single 
task Tcmb- 

3) Result: In the {N — l}-th iteration, the processor uti¬ 
lization waveform of and Tp is actually the processor 
utilization waveform of ti, ...,tm 
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As we can see, the scheduling of multiple tasks are treated 
recursively as the scheduling of two tasks. There are two 
challenges in the recursive algorithm. First, we need to model 
the processor utilization waveform as a single task Tcmb', 
second, we need to prove that the Assumption 13.11 and 13.21 
are satisfied during each iteration. 

First, we show how to characterize Tcmb as a single 
schedulable task from the processor utilization waveform. The 
processor utilization waveform within [ta{n),ta{n + 1)] have 
two modes. If 0 < z{n) < Ta{n), as illustrated in the upper 
figure of Fig. [3 the waveform within \ta{n)^ta{n + 1)] is 
divided into two segments. Each segment can be viewed as 
the execution of one instance of Tcmb- The task behaviors of 
the hrst instance is 

J Request Interval : ina.x{Ca(ji), z(ji)} 

Y Computing Time : C'a(n) + i?(n) 

and the task behaviors of the second instance is 

Request Interval : Tain) — max{C'tj(n), z{n)} 
Computing Time : (38) 

mm{Cj 3 {y{n)),Ta{n) - ma.x{Ca{n), z{n)}} 

On the other hand, if z{n) > Ta{n), as illustrated in the 
lower figure of Fig. [3 the waveforms within [ta{n),ta{n+l)] 
can be viewed as the execution of one instance of Tcmb, with 
the task behavior being 


A. Scheduling jnodel verification 

Consider the following scenario: 

1) Initially, two periodic tasks ri and T 2 are running on a 

single processor, with the following parameters: Ci = 
0.2min, Ti = Imin C 2 = 0.3min, Ti = 1.5min. ri 

starts at the time 0 while T 2 starts at the time O.Smin; 

2) An aperiodic task T 3 arrives at SOmin and another 

aperiodic task T 4 arrives at 50.6min. Both T 3 and T 4 
will stop after 57.Imin. Within [50.6,57.1], T 3 and T 4 
have the following characteristics T 3 (l) = 1 . 6 min, 
( 73 ( 1 ) = 0.5min, r3(2) = 2min, < 73 ( 2 ) = 0.6min, 
73 ( 3 ) = 1.7min, < 73 ( 3 ) = 0.2min, r3(4) = l.Smin, 

( 73 ( 4 ) = 0.4min, T 4 (l) = 2.5min, (74(1) = O.lmin, 

T4(2) = 3min, (74(2) = 0.4min, 24 ( 3 ) = Imin, 
( 74 ( 3 ) = 0.3min. 

3) At time llOmin, T 2 stops. Another periodic task T 5 ar¬ 
rives at 11 l.Smin, with T 5 = 1.2min and C 5 = 0.4min. 

4) A aperiodic task tq arrives at time llSmin and dis¬ 
appears after 120min. Within [113,120], tq has the 
following characteristics T6(l) = 4min, Tq{2) = 3min, 
( 75 ( 1 ) = 0.6min, Ce(2) = O.Smin. 

We are interested in the state of processor within [50, 57.1] 
and [110,120]. The waveforms shown in FiglHand Fig|3are 
consistent with those generated by Truetime. However, by 
using truetime, we have to initialize system whenever new 
tasks arrive and to simulate from the beginning. 


J Request Interval : Ta{n) 

\ Computing Time : Ca{n) + R{n) 

Next, we use mathematical induction method to prove the 
following Proposition. 

Proposition 6.1: In each iteration, for any time instant t, 
there exists fg and L such that the maximum request interval 
of Tq, is smaller than the minimal request interval of Tp within 

[to, to + L], 

Proof: In the first iteration, Tq denotes ri and Tp denotes 
T 2 . According to Assumption 13.11 and Assumption 13.21 the 
Proposition 16.11 holds. 

Suppose in the n-th iteration, for any time instant t, there 
exists fp and L such that the maximum request interval of Ta 
is smaller than the minimum request interval of Tp within 
[fo,fo + L], Tp denotes r„+i. Since Tcmb is derived from 
each request interval of r„, the maximum request interval of 
Tcmb is no larger than the maximum request interval of Tq. 
According to Assumption 13.11 we know that the maximum 
request interval of r„+i is smaller than the minimum request 
interval of r „+2 within [to, to + L], Therefore, the maximum 
request interval of Tcmb is smaller than the minimum request 
interval of t „+2 within [fo,fo + L] and Proposition 16.II in the 
{n + l}-th iteration holds. ■ 

VII. Simulations 

In this section, we verify the dynamic scheduling model by 
comparing task scheduling waveforms with those generated by 
TrueTime. Then, we simulate the scheduling and the battery 
models together to predict remaining battery capacity. 


Tasks Scheduling Within A Finite Time Window 



^^ 


50 51 52 53 54 55 56 57 

time(s) 

Processor Utilization Waveform Within A Finite Time Window 



time (s) 

Fig. 4. The State of Processor within [50, 57.1] 

As we can see in Fig|4]and|3 the value of a task will jump 
from 0 to 1 when it begins to execute. Then, the value of this 
task might go through several transitions between 1 and 0.5 
during the execution period. Finally, the value of this task will 
fall back to 0 once the execution is finished. 

B. SPSB simulation 

Fig. |6] shows the variation of battery capacity depending on 
processor status. When processor is busy, i.e. $(f) = 1, the 
battery capacity loss keeps increasing. When processor is free, 
i.e. <l>(f) = 0, the battery capacity loss begins to decrease due 
to the recovery effect. 

We simulate a battery same as that in ll22l . which has the 
parameters a = 40375 and = 0.273. We assume that the 
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Processor Utilization Waveform Within A Finite Time Window 


110 111 112 113 


115 116 117 11i 


Fig. 5. The State of Processor within [110,120] 


current drawn by the processor when it is busy is / = 200 
mA and the current vanishes when the processor is free. 

Fig. |6] shows the state of battery within two time windows 
[50,57.1] and [110,120]. According to Assumption 13.11 and 
13.21 we know that the state of the battery within its whole life 
period can be monitored by properly shifting the hnite time 
window. 



time (mintue) 


Fig. 6. Simulating the hybrid model within [50 57.1] 



Fig. 7. Simulating the hybrid model within [110 120] 


VIII. Conclusions 

In this paper, we have established analytical models for 
the behaviors of multiple aperiodic tasks scheduled on a 
single processor supported by a single battery. We assume 
that the tasks are scheduled under a RMS like algorithm 
that assigns priority of tasks based on information within a 
hnite time window. Our model is presented using a set of 


nonlinear difference equations that describe the task behaviors 
together with continuous differential equations that describe 
battery discharge behavior. One may hnd that these equations 
can be studied as a hybrid system. It can be perceived that 
through our modeling effort, battery behavior and scheduled 
task behaviors can now be studied jointly within the same 
theoretical framework. This fact is well aligned with the goals 
of cyber physical systems theory. 
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